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In this paper, a purification protocol is presented and its performance is proven to be optimal 
when applied to a particular subset of graph states that are subject to local Z-noise. Such mixed 
states can be produced by bringing a system into thermal equilibrium, when it is described by a 
Hamiltonian which has a particular graph state as its unique ground state. From this protocol, 
we derive the exact value of the critical temperature Tcrit above which purification is impossible, 
as well as the related optimal purification rates. A possible simulation of graph Hamiltonians is 
proposed, which requires only bipartite interactions and local magnetic fields, enabling the tuning 
of the system temperature. 



Introduction. Any quantum technological implementa- 
tion is plagued by environmental noise. The possibility 
to purify quantum states, or to use error correcting al- 
gorithms to stabilize quantum operations, is therefore a 
necessary step towards reaping the benefits of quantum 
technologies. Much attention has been focused lately on 
the purification 0, 0i 0i S S IE of a large class of mul- 
tipartite states called graph states 0. After the initial 
restriction to two-colourable graph states J 0| , the ideas 
have been extended to all graph states [7|, and general- 
ized to other stabilizer states . The variety of differ- 
entprotocols trade off between a large tolerance to noise 
0, |i| and the rate of purification . Graph states have 
proven important for realizing a variety of quantum infor- 
mation tasks such as performing quantum computation 
jlOl | , quantum communication and as a means for effi- 
ciently approximating other quantum states, such as the 
ground states of strongly correlated systems ■ While 
some bounds have previousl y be en proven on the ability 
to purify multipartite stat es Il,'ill4l , optimal results only 
exist for two-qubit states llsLllq . 

Here we concentrate on the purification of graph states 
that are subject to the physically motivated independent 
Z-noise. From the technological perspective, large-scale 
quantum computation is still too difficult to implement. 
However, graph states can be made and manipulated in 
the laboratory, e.g. by controlled collisions of alkali atoms 
trapped in optical lattices 0i Q ■ To enable the con- 
trolled collisions, it is necessary to employ two magnet- 
ically sensitive hyperfine levels of the atoms to form a 
qubit. This sensitivity means that the states are subject 
to decoherence from stray magnetic fields, in addition to 
any uncontrolled collisions that may occur. These errors 
are described by Z-errors, and as experiments improve, 
one can expect them to become more localised. An al- 
ternative approach to preparing a graph state involves 
implementing a Hamiltonian, known as a graph Hamil- 



tonian, which has the desired state as its ground state. 
As it is impossible to cool the system to absolute zero, 
the resulting equilibrium state will always be a thermal 
state. For graph Hamiltonians, the thermal noise corre- 
sponds to local Z- noise on each qubit. 

In this paper, we consider a certain purification proto- 
col applied to arbitrary graph states in the presence of 
independent Z-errors. Whilst this multipartite protocol 
is not novel or sophisticated 0, lisl , it has the advan- 
tage of being analytically tractable. Most importantly, 
we prove optimality of this protocol for these types of 
error, both in the sense of the level of noise that can be 
tolerated as well as the scaling of the purification rate, 
for a specific subset of states that includes the cluster 
and Greenberger Horne Zeilinger (GHZ) states. For up 
to seven qubits, this subset can be shown to be isomor- 
phic to arbitrary graphs under local operations. We pro- 
pose a method for simulating the graph Hamiltonians 
that is comprised of only two-body collisions and local 
magnetic fields, even though the resulting Hamiltonian 
has at least 3-qubit interactions |2oLl2]|. The critical pu- 
rification temperature related to these models can easily 
be made to lie above the typical temperatures given, e.g. 
from optical lattice realizations of such Hamiltonians. 

Graph Hamiltonians and Graph States. Let us intro- 
duce a graph G, which is defined by a set of vertices 
Vg, and a set of edges Eg, describing the connections 
between the vertices. To each vertex of this graph, we 
attach a spin-i particle (qubit), and define a graph state 
to be the ground-state of the following Hamiltonian, 
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The Si's are the coupling strengths, which we henceforth 
take to have equal magnitude, Bi ~ B, and we assume 
B > 0. The interaction terms Ki commute with each 
other, [Ki,Kj] = 0, and hence each term individually 
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FIG. 1: Purification protocol for independent noise. First, 
we take many copies of the noisy graph state, (a) We con- 
struct from these two-qubit nearest-neighbour states (noisy), 
(b) Two-qubit states are purified (if possible) . (c) ControUed- 
phase gates are applied between local qubits. (d) All qubits 
except one from each party are measured in the X-basis, leav- 
ing the remaining qubits in the purified state. 

stabilizes the eigenstates oi H. A local Pauli Z- rotation 
on qubit i, Zi, commutes with all the Kj where j ^ 
and anti-commutes with Ki. Hence, the excitations of 
the Hamiltonian are given by local Z-rotations applied 
to the ground state. 

A constructive way to produce graph states is found 
by close analogy with the cluster states |^ . We can 
produce a general graph state by creating the |+) = 
(|0) -I- |1))/a/2 state on each vertex of G, and applying 
controUed-phase gates along the edges of the graph. The 
action of measurements on the cluster state also applies 
to the graph states, such that when a Z-measurement is 
performed, all nearest-neighbour bonds are severed. This 
means that graph states can be cut into sections using 
Z-measurements, and can be combined together using 
controUed-phase gates. The key to what follows is the 
realization that Z-errors commute with both of these op- 
erations, and consequently remain as Z-errors. 

Examples of graph states that are of particular inter- 
est in quantum computation and communication are the 
cluster and GHZ states. The first state corre- 
sponds to a graph that is given by a d-dimensional cube, 
while the latter state corresponds to a graph that is lo- 
cally equivalent to a single qubit connected with all other 
qubits. 

Purification of Thermal Graph States. The thermal 
state of the Hamiltonian in Eqn. of |Vg| — N qubits 
at temperature T = l/{kBP) is given by 
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where we set the Boltzmann constant equal to unity. 
This density matrix can be written in terms of local Z- 
errors as p{p) — £i£2--£n IV') (V'L where 

£iP = [(1 - p)p + pZipZi\ 
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is the probability of a Z-error occurring at a certain site 
due to the non-zero temperature T. The graph state |V') 
is the unique ground state of Hamiltonian . Our aim is 
to purify towards the state \ip) using many copies of p{p)- 
We consider that each vertex of the graph is controlled 
by a different party and that operations such as mea- 
surements and controUed-phase gates are only allowed 
locally, but involving the many copies. This restriction 
corresponds to the scenario of quantum repeaters [23| . 
where the different parties are physically separated, and 
also serves to illustrate the entanglement properties of 
the system. We do not envisage implementation of the 
protocol in other scenarios, since it is generally cheaper 
in terms of resources to create the state directly - our 
aim is to prove optimality, providing a benchmark for all 
other protocols. 

The purification procedure consists of breaking down 
the graph states into smaller blocks, purifying them, and 
then recombining them. We refer to it as the Divide 
and Rebuild Purification Protocol (DRPP). The smaller 
blocks that we choose to use are two-qubit states, for 
which there are analytic purification results, returning 
maximally entangled states. As already specified, the 
splitting of the graph state is readily achieved with Z- 
measurements. This leaves us with a two-qubit mixed 
state P2ip), which can be purified to a maximally entan- 
gled state 1-02) provided 

Once we have generated a maximally entangled state 
for each of the edges of the graph, each local party (i) 
holds a number of qubits equal to the number of nearest- 
neighbours \Eq\- These can be reduced to a single qubit 
that contains all these links by applying controUed-phase 
gates and performing AT-measurements, as described in 
Fig. [2 

The condition for purification of the two-qubit state 
between neighbouring sites is readily found to be 
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This is known to be necessary and sufficient for purifying 
states that are Bell diagonal, like the states considered 
here 0, 0, 01 ■ We can hence calculate that the maxi- 
mum temperature at which the DRPP can purify graph 
states is given by 
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Optimality and Rates. While the DRPP can be ap- 
plied to any graph state (with minor modifications), we 
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FIG. 2: If we assume the multipartite state can be purified, 
then this implies that we can purify the two-qubit state. Con- 
versely, if the two-qubit state cannot be purified, the assump- 
tion must be broken, (a) Alice and Bob take the two-qubit 
state and reconstruct the noisy graph state, (b) This state 
is purified by assumption, (c) All extra qubits are measured 
out to return the original pair, now pure. 



can prove its optimality for a specific subset of graphs. 
For clarity, we shall restrict to only cluster states (of ar- 
bitrary dimension) and GHZ states, all of which fall into 
this classification. The optimality that we prove is with 
respect to both the level of noise that can be tolerated 
such that purification is still possible, and with respect to 
the number of copies of the initial, noisy, state required to 
form a single pure copy. In particular, one can prove that 
states with higher levels of noise than the one dictated by 
Eqn. Q can never be purified with any protocol. We do 
this by considering a purification protocol for a two-qubit 
state. The two parties sharing this state are allowed to 
introduce extra qubits and the local operations that they 
apply subsume multipartite considerations for the addi- 
tional qubits. 

If two parties, Alice and Bob, hold several copies of the 
noisy two-qubit state P2{p) .25], they can locally recreate 
the initial thermal state, p{p)- For linear graphs (Fig.|2J), 
this simply corresponds to Alice and Bob locally creating 
their own thermal cluster states, and connecting them 
to p2{p) with controlled-phase gates. In the case of a 
more general graph, we require some additional connec- 
tions between Alice and Bob. These are achieved by us- 
ing multiple copies of P2{p) (Fig. E^b)). At this stage, 
we assume that purification of p{p) is possible, yielding 
From there, Alice and Bob can measure out all the 
qubits that they added, leaving a pure two-qubit state. 
Hence, if p{p) can be purified, P2{p) can always be puri- 
fied. However, we know that p2{p) cannot be purified if 
(1 — pY < ^. Hence, under this condition, our assump- 
tion must be false i.e. the multipartite state p(p) cannot 
be purified, whatever the protocol. The DRPP saturates 
this bound and hence is optimal for independent Z-noise. 
While analysis of the protocol of 0, is a difficult prob- 
lem, numerical results indicate that it also saturates the 
bound, and hence is optimal for this type of noise. 

The rate of purification, i?^, of the graph state 1-0) can 
be calculated in terms of the yield of a Bell state, i?2- We 
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FIG. 3: (a) Simplest example of a 2ci-cluster state, shared be- 
tween 4 parties, (b) Alice and Bob can locally reconstruct the 
square graph using two copies of p2, and applying controlled- 
phase gates between them, (c) For the triangular configura- 
tion, the optimality proof fails. 
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FIG. 4: By measuring in the Z basis on every third qubit, 
we can create two-qubit states between all nearest-neighbours 
with only three copies of the original linear cluster state. 
Hence, we say A'gco = 3. 

take the standard definition of rate. 
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If we can purify p into at a rate (i.e. we require 

copies of p to create IV')), then we can create a 
Bell state between any linked pair just by performing Z- 
measurements on excess qubits. There could be a more 
efficient way to generate this Bell pair, requiring fewer 
copies, so i?^ < i?2- 

Similarly, if we can purify Bell pairs, then we can gen- 
erate For that we take a Bell pair between each 
nearest-neighbour and the local parties perform the re- 
construction as specified by the DRPP (e.g. controlled- 
phase between local states and AT-mcasurements). Most 
of the Bell pairs can be purified in parallel - we only re- 
quire iVgeo copies of p to purify enough copies, where the 
geometric factor TVgco depends only on the local degree 
of the graph, \Eq\: and is otherwise independent of the 
number of qubits in the graph. This is because the Z- 
measurements commute with the errors, and divide the 
state p into separate blocks. As a result, R^ > 



Combining the two results. 



i?2 > i?V > 



gco 



For c?-dimensional cluster states, we can readily evalu- 
ate iVgoo , since we must measure out all qubits connected 
directly to the pair that we are interested in isolating. 
Starting from the edge of the lattice, this uses up 3d — 2 
qubits (we include in this number some single qubits that, 
while we might not intend to measure them, become iso- 
lated), plus the two for the state p2- If our lattice extends 
to N qubits in each direction, then we can generate on 
average (A^ — l)iV''^^/(3o?) copies of p2 from a single copy 
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of p. We need d{N — l)N''-~^ different copies of p2, and 
hence iVgoo = id? ■ As demonstrated in Fig. 01 this cor- 
responds to A'goo = 3 for d = 1. The resulting rate is 
independent of the number of qubits in the system, only 
depending on R2, which is optimal up to a small numer- 
ical factor. For A^-qubit GHZ states, A^gco = — 1. 

Physical Implementation. A significant achievement 
in recent quantum engineering experiments is the con- 
struction of cluster states with optical lattices They 
are produced with a single operational step independent 
of the size of the system. It is natural to consider this 
setup for studying the purification properties of thermal 
cluster states. This requires the implementation of both 
the cluster Hamiltonian and a purification protocol in a 
physical setup. In the following, we will present a simple 
way of simulating the Hamiltonian for cluster states, 
using proven experimental techniques in optical lattices. 

When an entangled state decoheres, there is a charac- 
teristic lifetime that determines when the state becomes 
separable. On the other hand, if it is possible to ob- 
tain an interaction described by the Hamiltonian that 
has this entangled state as a ground state, then provided 
the energy gap is large enough in comparison to the de- 
coherence rate, entanglement can survive indefinitely in 
the system, e.g. in the form of purifiable mixed states, as 
we have already seen. 

Simulating Hamiltonian for a general graph is a rel- 
atively straightforward task. The method we adopt here 
consists of a unitary operation, Uq \Ug\ , applied before 
[after] the evolution with respect to a local Hamiltonian. 
This evolution is generated by applying a uniform mag- 
netic field in the x-direction. When Ug corresponds to 
controlled phase gates between all pairs of qubits con- 
nected in the graph, then it is easy to show that the 
resulting effective Hamiltonian 

H = UG(^Bj2x}jul, (6) 

is of the form |^ [2^ . The spectrum of the Hamiltonian 
corresponding to the magnetic field B ^ ■ Xi is the same 
as that of the Hamiltonian H, as they are related by an 
isospectral transformation. Hence, the generated thermal 
state directly corresponds to the one of H. 

In optical lattices, the unitary Ug can be realized in 
cubic lattices by controlled collisions between nearest- 
neighbours. This is preci sely the operation which is ex- 
perimentally employed in [l7| for the generation of a clus- 
ter state. Meanwhile, the local magnetic fields are imple- 
mented by globally applied Raman transitions between 
the hyperfine states that encodes a qubit. Thus, the re- 
alization of the graph Hamiltonian is readily achieved. 
The stationary state of this system is the thermal state 
of Eq. 0. Further, this system gives us the abi lity 
to vary the temperature. In previous experiments |29| 
the recorded temperature after performing optical cool- 



ing was given by the relation T « 0.1 Co, where Uq is 
light shift potential created by the optical lattice. From 
the value of the critical temperature given in ||SJ|, one de- 
duces that by employing moderate local magnetic fields 
with amplitude B > O.lC/o, one can bring the system into 
the purifiable regime. 

Once we are able to implement a purification proto- 
col in optical lattices, whether it be the DRPP or any 
other, this provides us with the perfect test-bed to probe 
the maximum temperature that still allows purification, 
and verify the critical temperature given in Eqn. [S] The 
potential implementation of such schemes has been de- 
scribed extensively in, for example, j^. The main draw- 
back is the requirement of local addressability, which can 
be circumvented with the help of superlattices [ll . E^E^ 
by breaking of the translational invariance of the lattice. 

Conclusions. Here, a purification protocol has been 
proven to be optimal when applied to a sub-class of graph 
states, including the cluster and GHZ states, subject 
to Z-errors. Although we have restricted ourselves to 
this particular form of noise in this paper, we emphasize 
that thermal states arising from graph Hamiltonians are 
precisely of this form, rendering the considered types of 
states interesting also from an experimental and practical 
point of view. 

While the optimality proof cannot be applied to cer- 
tain types of graphs (e.g. the icosahedron) , and we are 
aware of examples of noise (e.g. local or global white 
noise) where the proposed protocol is not optimal, sev- 
eral extensions of our results are possible. Two of the 
authors have shown [^O] that the optimality proof can 
be extended to a wider class of states and to different 
forms of noise. For instance, it can be shown that all 
graph states of up to seven qubits can be brought by 
local unitary operations to a form where the optimality 
proof can be applied 0- Thus, all these graphs have a 
critical temperature given by Eqn. In addition, it 
shown in (20] how the optimality proof presented in this 
paper can be extended to other types of noise, providing 
an upper-bound to the error probability that can be pu- 
rified. Further extensions taking into account non-graph 
states and noisy local operations during the purification 
protocol will also be examined. 
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